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Grand unified theories may display multiply interacting fields with strong coupling dynamics. 
This poses two new problems: (1) what is the nature of chaotic reheating after inflation, and (2) how 
is reheating sensitive to the mass spectrum of these theories ? We answer these questions in two in- 
teresting limiting cases and demonstrate an increased efficiency of reheating which strongly enhances 
non-thermal topological defect formation, including monopoles and domain walls. Nevertheless, the 
large fluctuations may resolve this monopole problem via a modified Dvali-Liu-Vachaspati mecha- 
■ nism in which non-thermal destabilsation of discrete symmetries occurs at reheating. 
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An ideal inflationary scenario should arise naturally from within supergravity or a Grand Unified Theory (GUT) 
without fine-tuning. It should solve the plethora of problems of standard cosmology while simultaneously diltuting 
the monopoles inevitably produced due to the nomotopic content 7r 2 (G/?7(l)) ~ 7Ti(J7(l)) ~ Z of the standard 
ly-j , model. Significant progress has been made within GUT's and supersymmetric theories towards this utopic vision 
0,^). However, the issue of reheating after inflation in these theories, where the universe is revived, phoenix-like, from 
the frozen vacuum state, has remained relatively unexplored || . This is precisely one area where the full symmetry 
and particle content of the underlying theory is likely to be crucial. 

Indeed reheating poses a severe threat to the simple ideal presented above since non-perturbative effects are typically 
dominant 0-0]. Reheating is therefore not a minor phase at the end of inflation, of little dynamical interest. The 
large quantum fluctuations allow for GUT baryogenesis and, as pointed out by Kofman et al B, may cause 
the monopole and domain wall problems to reappear due to non-thermal symmetry restoration. This last possibility 
is actually rather difficult in simple models of reheating with only two fields Q. As we shall show, however, this 
situation changes dramatically in the case of multiple fields, relevant for GUT models. 

The main motivation of this work then is to understand what new effects multiple fields have on reheating. This 
issue encompasses two particularly interesting unknowns, (i) The nature of reheating at strong coupling when the 
fields evolve chaotically. This is relevant for GUT's with divergent UV fixed points [p|JlO|] and models such as softly 
broken Seiberg-Witten inflation ]TT[ | where reheating occurs in the strongly coupled, confining, regime. Setting aside 
the subtle issue of the quantum behaviour of gauge theories at strong coupling ]l2j , in the two and three scalar-field 
cases studied so far, (classical) chaotic motion has been typical jl3|], especially at reheating. This chaotic evolution 
parallels results in the Einstein- Yang-Mills equations [O], semi-classical QCD and lattice gauge theory Jlj|. Thus a 
natural question is "what is the nature of chaotic reheating ?" 

The second issue is (ii) the sensitivity of reheating to the mass spectrum of the theory. For example, in SO(10), the 
one-loop effective potential is built from four quadratic invariants |l6[ ]. "How does reheating depend on the relative 
masses of these four fields ?" 

We find that both these factors can significantly enhance the efficiency and power of reheating. This in turn im- 
proves the possibilities for GUT baryogenesis and non-thermal symmetry restoration (NTSR), with the concommitant 
formation of topological defects |J] . 

Consider the inflaton condensate <f> in the reheating phase with large couplings A* to the fields <5 . Further, couple 
<j> to the minimally-coupled scalar field \ p/fl] via the interaction term \g 2 <fi 2 X 2 1 where g is a dimensionless coupling 
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constant. Then the modes of \ obey |10[ : 



d 2 (a 3 / 2 Xk) , fk 2 2 3k _ 2 2 ., ., 

2+ m x + -rP + s^ )( a Xk) = o (i) 



dt 2 v« 2 x 4 

where a(t) is the scale factor of the universe obeying the Hamiltonian constraint H 2 = (a/a) 2 = ptot/3 with p to t the 
total energy density and p the total pressure of the universe, k = 8irG = 8tt. M p i is the Planck mass. 

In the simple model where the *S? fields are absent, p — <j> 2 /2 — V(<fi), and if the inflaton effective potential is 
V(4>) = to|0 2 /2, the solution is || (j) — ®(t) sinm^t, where <& decays slowly with the mean expansion. 

In this case eq. (|l|) can be cast in the form of the Mathieu equation 0: (a 3 ^ 2 Xk)" + [A + F](a 3 / 2 x,k) = 
with F = — 2qcos(2t'). The adiabatically evolving parameters Q A — k /(rrv^a 2 ) + m 2 ,/m'^ + 2g 2 $ 2 /m 2 j and 
q = $ 2 (g 2 /rr& + 3k/16), span a plane geometrically dissected into stability and instability regions (fig. lb) ||. The 
solutions to the Mathieu equation in the instability regions are periodic with envelope y ~ exp(/xm0i), and Floquet 
index \i > 0. This resonance continues until backreation due to ((^x) 2 ) shuts off the resonance by forcing all physical 
(k > 0) modes out of the dominant first resonance band 

Now switch on the couplings A* to the fields \I> and move into the strongly coupled, chaotic region of the parameter 
space. While we are unable to study chaotic reheating in full generality, we have full control over the region in which 
the chaotic fluctuations are extremely rapid. This is because of the following little-known theorem (JT8): In the limit 
of rapid variations, chaotic flows become uniformly indistinguishable from white noise Jl9| |. This is an extremely 
powerful result since the investigation of reheating becomes very pure and tractable. 



In this limit, eq. (|l|) again reduces to quasi-Mathieu form, but now with a new potential F |20 21||: 

F = -2qcos(2t)+g 2 £(t). (2) 

Eq. corresponds to A*-independent stochastic evolution Gaussian white noise) describing the effect of 

multiple fields in the strong coupling limit (in which case q — 0, as there is no periodic component). This form of the 
potential also models the backreaction of quantum fluctuations on the mean periodic evolution of <$> (q ^ 0) in the 
case where the \t are absent |p2| , po| . 

Using spectral theory it was shown that the stability bands essentially disappear in the case of stochastic 
potentials F. In the strong-coupling limit, g — > oo we have a rather intimidating analytical result for /i^ which 
is rather intractable and strongly emphasises the need for numerical study of eq. (|l|) . A selection from the extensive 
simulations performed [j23f for this problem are shown in figs. (2a, b). The main generic results are the complete 
breakup of the stability bands of the Mathieu equations, with \i > Vfc, and the significant increase of fi, over the 
purely periodic case. This is seen in all our simulations and is illustrated in figs. (2a) and (3a) which give fi along the 
physical separatrix (k = 0): A ~ 2q. 

A crucial result of these large fi is the possibility of non-thermal symmetry restoration (NTSR) [Q. Consider the 
simplest 2-ficld effective potential with symmetry breaking: 
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V(^ X )=--{^-cfl) 2 + y -cf X 2 . (3) 

This gains the following quantum corrections: AV = 3\((5(Q 2 )<f> 2 /2 + g 2 ((Sx) 2 )4> 2 /2 + <7 2 (((5c/>) 2 )x 2 / 2 B> where the 
dominant variance is given in the Hartree approximation by M {(S(f>) 2 ) = (2n 2 a 3 )^ 1 j dkk 2 \a 3 / 2 5(f)k\ 2 ~ / <ifcfc 2 e 4 ' im * t . 
From energy conservation, one gets the bound Jl {(5(f>) 2 ) < Cg" 1 XM p i ln~ 2 g~ 2 with C ~ 10 -2 . Saturating this bound 
implies that Eq. (J3J) can have a positive effective mass, e ^ = (V + AV)" even for 0o ~ W 16 GeV. 

Nevertheless, exhausting the available energy and realising these large variances in the simple Mathieu equations 
is highly non-trivial jf|: (i) sufficiently large [i are required since {(5(f) 2 ) ~ J dkk 2 e 2tlt . (ii) The breadth (and very 
existence of) the first instability band controls how long the resonance continues before backreaction shuts it off (see 
fig lb), (iii) The expansion of the universe is forced to be monotonically decreasing: H = —n(j) 2 /2, which drives 
(A, q) — > (0,0), damping the resonance and reducing the variances |p4f . (iv) Finally, for defect production we require 
e j:j > for time scales St > uj^ 1 , the period of <fi oscillations. 

These factors conspire against NTSR and defect production in the simplest models of preheating However, in 
the strong-coupling limit of GUT's considered here we know that (i) we can easily achieve very large /x (see fig. 2a), (ii) 



since the stability bands are completely destroyed |20 21]] , modes never stop growing, fi > 0, even when backreaction 
becomes important, (iii) Finally, and very importantly, when there are multiple fields, the expansion need not be 
monotonic and indeed, in our case will increase and decrease stochastically ]l3j until reheating and thermalisation are 
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completed. This means that the expansion need not only have a damping effect, but can also act as a pump, when 
H > 13 1 , enhancing the fluctuations ((8cj)) 2 ) and ((<5x) 2 )- 



From these robust considerations, we expect NTSR to be much more effective at strong-coupling, allowing us 
to approach the limiting variances set by energy conservation. Indeed this is bourne out by direct simulation of 
m 2 eff cx \{(8<j>) 2 ) + g 2 {(8 X ) 2 } (see figjj). 

Examining the statistics of m? e ^ |23| we see that while e ^ fluctuates very rapidly, very large values of ((<5x) 2 ) 
occur, with rrA e tf S [— A0q, g 2 M 2 t ). Assuming \(j) 2 -C g 2 M 2 t , the statistics of rrA e ff are highly skewed, yielding 
rn^ > 0. This means that most of the time symmetry is restored and <f> is likely to diffuse across the origin, leading 
to defect production including monopoles in a full model including 5C/(3) C x SU(2)l x U(1)y- 

We now study the effect of mass spectrum deformations on reheating ||| . A simplified toy model to study this 
issue is given by the 3-field effective potential: V{<p,ip,x) — m 2 cf) 2 /2 + m 2 ip 2 /2 + g 2 (f) 2 x 2 /2 + A 2 </? 2 x 2 /2. Eq. (|l|) is 

F gp = —q(cos 2t + cos 2nt) . (4) 



modified and leads to a new F |21LE5 



where we have chosen m^/m^ = tt, which implies that (|4|) is a quasi-periodic function []21 23 since the masses are 
irrationally related. 

In our case, spectral theory results guarrantee that the stable bands generically form a Cantor set ^ij, which are 
often of very small measure - hence physically unimportant - mimicking the stochastic case presented earlier. Again 
however, no estimates for n are available analytically. Our numerical results show the growth of n in this case and 
the significant widening and steepening of the instability bands relative to the Mathieu case (compare figs. 1 and 4). 
Also interesting is the development of "non-thermal edges" in fi (fig. 4a and p3fl). 

From our discussion of the stochastic case, we expect ((S(f>) 2 ) to increase and again this is bourne out numerically. In 
figs. (5a, b) we plot log((<50) 2 } for the pure Mathieu and quasi-periodic cases as functions of q. The maximum variance 
is ~ 1.5 orders of magnitude larger in the quasi-periodic case after the short time t = 23m < ^" 1 . The quasi-periodic 
model therefore is significantly more efficient at restoring symmetry. In general the NTSR strength of GUT theories 
will depend sensitively on the mass spectrum of the theory, here encoded by m^/m^. 

The monopole problem has now again become a major concern in large ( "chaotic or incommensurate" ) regions of 
the coupling/mass parameter space, including the expansion of the universe. Due to the non-thermal, quench-like, 
nature of the symmetry breaking, the correlation length 5 of the fields will be much smaller than in the equilibrium 
case, and therefore the defect density oc S~™ in — 1 for domain walls, n = 3 for monopoles) will be correspondingly 
larger than the equlibrium Kibble prediction J26| . 

If NTSR succeeds, a second stage of inflation will occur j|,^7j while the vacuum energy V(0) dominates over 
the energy of the ((Sx) 2 )- During this time a(t) increases by a factor ~ (cP/A) 1 / 4 for the potential (||), which 
cannot, therefore, supply the needed ~ 20 e-foldings to dilute the monopole density sufficiently In 50(10) or 
5£/(6), however, the monopole transition is separated from the lower transitions and this may allow enough secondary 
inflation to dilute the monopole abundance sufficiently. 

However, this is rather model dependent and the large corrections to the effective potential offer us an alternative 
escape route via defect-defect interactions. The full corrections to the effective potential, include not only the quadratic 
contributions affecting m| e j f , but also odd powers which do not respect any previously existing discrete symmetries 
(required for example for successful D-term inflation p8[ ). The odd-power corrections to eq. (^), with g = 0, are: 

AV odd = \<j)(5f - (j) 2 5(j)) + A<5# 3 (5) 

These terms softly break the Z2 4> — > —<j) symmetry of eq. (|J). This allows an implementation of the Dvali-Liu- 
Vachaspati mechanism for solving the monopole problem as follows: imagine that NTSR was successful enough to 
produce monopoles and domain walls during reheating. The 8(f) terms in eq. (|5|) automatically cause the domain walls 
to be unstable p9[ | since the minima at ±<fio are no longer degenerate, causing a pressure difference ~ 2AV dd(4>o) 
across the walls. 

Since after preheating <fi <C ({8(f>) 2 ), the first term of (||) dominates. The constraint that the walls percolate gives 
us (8(f)) < 10 _1 A/ p ;. Monopoles are then swept up on the walls and dissipate because the full symmetry (e.g. 5t/(5)) 
is restored there |30|. Requiring that the pressure difference drives the domain walls to collapse and decay before 
dominating the energy density of the universe yields @ (8<p) > W^X^^M^ - l()- 3 M p i if A - 10~ 3 . These 
bounds are exactly in the range of values expected if NTSR is successful. The advantage of this formulation is that it 
provides a specific implementation of the general Dvali-Liu-Vachaspati mechanism and uses the same large quantum 
fluctuations which produced the monopoles to remove them. 
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FIG. 2. Stochastic eq. (2), g = 5: (a) Floquet index /i vs q along A — 2q. (b) Contour plot of instability chart - note the 
breakup of Mathieu bands and the very large peaks in /i. 
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FIG. 3. Three realisations of ((<5x) 2 } vs g at t = 2m^ 1 for q = 5. At larger times, {(8x) 2 } becomes completely dominated by 
a single peak (c.f. fig. 2(b)). 
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FIG. 4. Quasi-periodic case, eq: (4) (a) /i vs q on A — 2q. (b) Instability chart (A, q). Note the larger n and proliferation of 
instability bands compared with figs. (1 a,b). 




FIG. 5. log((5(j>) 2 } vs q for (a) the Mathieu eq., and (b) the quasi-periodic eq., both calculated at t = 23m^ > 1 . The maximum 
value of {(8(f)) 2 } is about 20 times larger in (b). 
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